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ABSTRACT 

Context. The gamma-ray binary LS 5039 and the isolated pulsar PSR J1825-1446 were proposed to have been formed in the super- 
nova remnant (SNR) G016.8-01.1. 

Aims. We aim to obtain the Galactic trajectory of LS 5039 and PSR J1825- 1446 to find their origin in the Galaxy, and in particular 
to check their association with SNR G016. 8-01.1 to restrict their age. 

Methods. By means of radio and optical observations we obtained the proper motion and the space velocity of the sources. 
Results. The proper motion of PSR J1825-1446 corresponds to a transverse space velocity of 690 km s -1 at a distance of 5 kpc. 
Its Galactic velocity at different distances is not compatible with the expected Galactic rotation. The velocity and characteristic age 
of PSR J1825-1446 make it incompatible with SNR G016. 8-01.1. There are no clear OB associations or SNRs crossing the past 
trajectory of PSR J1825- 1446. We estimate the age of the pulsar to be 80-245 kyr, which is compatible with its characteristic age. 
The proper motion of LS 5039 is fi a cos 6 = 7.09 and ji s = -8.82 mas yr" 1 . The association of LS 5039 with SNR G016. 8-01.1 is 
unlikely, although we cannot to discard it. The system would have had to be formed in the association Ser OB2 (at 2.0 kpc) if the age 
of the system is 1.0-1.2 Myr, or in the association Set OB3 (distance 1.5-2 kpc) for an age of 0.1-0.2 Myr. If the system were not 
formed close to Ser OB2, the pseudo-synchronization of the orbit would be unlikely. 

Conclusions. PSR J 1825- 1446 is a high- velocity isolated pulsar ejected from the Galaxy. The distance to LS 5039, which needs to 
be constrained by future astrometric missions such as Gaia, is a key parameter for restricting its origin and age. 

Key words, stars: individual: LS 5039 - pulsars: individual: PSR J1825-1446 - radio continuum: stars - proper motions - X-rays: 
binaries - gamma rays: stars 



1. Introduction 

During a core-collapse supernova explosion of a massive star, 
the remaining compact object, either isolated or in a binary sys- 
tem, can receive a natal kick and obtain a h igh peculiar veloc- 
ity (iHillsl 11983b I van den Heuvel et al.l l2000h . A binary system 
can obtain a m oderate velocity from a symmetric Blaauw kick 
(Blaauw 1961), whereas the compact object can acquire addi- 
tional momentum from an asymmetric supernova kick dStonel 
1982). If the system remains bound after the explosion, the 
parameters of the binary system before and after the super- 
nova explosion can be related by measuring the high pecu- 
liar velocity of the system, in parti cular for X-ray binaries (see 
e.g. Brandt & Podsiadlowski 1995, and the particular cases in 
iTauris et al.l Il999b iMartin et al.l 12009b . On the other hand, the 
space velocity of a binary system can be used to determine its 
past Galactic trajectory (see e.g. iRibo et ai1l2002t iDhawan et alJ 
120071 iMiller- Jones et al.ll2009ah . Of particular interest are the 
runawa y systems, which have acq uired very high peculiar ve- 
locities dMirabe l et al. 2001, 2002). With enough information, a 
complete evolutionary history of the binary system can be ob- 
tained with constraints on th e progenitor and the formation of 
the co mpact object, (see e.g. IWillems et alJl2005l : iFragos et alJ 
2009). Among the binary systems receiving natal kicks we can 
find the special case of ga mma-ray binaries, wh ich usually form 
eccentric binary systems (ICasares et al.ll2012bl) . 



Only a few binary systems have been associated with 
gamma-ray sources, and therefore it is important to understand 
how they formed, and how they are distributed in the Galaxy. 
However, the origin of most of the binaries w ith gamma-ray 
emission is still unknown. Mirabel et al .1 (12004 studied the ori- 
gin of the gamma-ray binary LS I +61 303, and IDhawan et al.l 
(2006) measured a precise proper motion that allowed the au- 
thors to obtain the space velocity of the source and constrain the 
mass lost in the supernova event that formed the compact object 
of the system to ~1 M . The natal kick of the supernova that 
produced th e pulsar in the sys tem PSR B 1259-63/LS 2883 was 
discussed in lHughes & Bailesl i 19991) and lWang et ail d2006l) . al- 
though no proper motion has been measured for this system. 
The Galac t ic m otion of LS 5039 was ex tensively discussed in 
IRibo et al.l d2002l) (see details in Sect. lQilMirabel & Rodriguesl 
(120031) and lReid et alJ d201 ll) used high-accuracy radio astro m- 
etry to compute the Galactic velocity of Cygnus X- 1 , which al- 
lowed them to obtain the mass of the progenitor of the black 
hole, and its association with the Cyg OB2 association. They 
also determined that the black hole was formed without a su- 
pernova explosion. For Cygnus X-3, IMiller- Jones et al. (2009b) 
obtained an accurate proper motion of the source that was used 
to better identify the position of the X-ra y binary core, although 
the origin of the system is still unknown. Mold on et al.l (|201 lb) 
determined for HESS J0632+057 that its total proper motion is 
below 4 mas yr _1 . Other recent gamma-ray binary candidates 
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are 1FGL J1018.6-585 6 (ICorbetet alJ20lU iPavlov et al.ll2011b 
Ackermann et al. 2012) and AGL J224 1+4454 (I Williams etaD 
2010; lCasares et al.l2012bl) . However, the propermotion of these 
sources is still unknown. These examples show the potential of 
a good determination of the proper motion of a system, and also 
the need for better astrometry in this peculiar binary population. 

On the other hand, accurate astrometry of isolated pulsars 
opens a wide field of scientific research including fundamen- 
tal reference frame ties, physics in the core-collapse supernovae 
and imparted momentum kicks, association with supernova rem- 
nants (SNRs) and determination of the age of pulsars, model- 
independent estimates of distances through parallax measure- 
ments, or determination of the distributio n of electron density in 
the interstellar med ium (see for example Ichatteriee et al.l l2009: 
Brisk en et al.l 12002 and references therein). It is common that 
pulsars acquire high velocities at birth, becoming the fastest 
population in the Galaxy, with a mean t ransverse velocity of ~ 
450 km s" 1 , and up to above 10 3 km s" 1 (Lvne & Lorime rll994t 
iHobbs et al.l l2005t IChatterjee et all 120051) . The pulsar velocity 



distribution provides information abou t the supernova sym me- 
try and the binary population synthesis (Bris ken et al.i l2003). 

Here we present the results from two high-resolution astro- 
metric projects to determine the proper motion of LS 5039, and 
the nearby pulsar PSR J 1825 -1446, the only known compact 
objects in the field of the SNR G016.8-01.1. In Sect.|2]we de- 
scribe the gamma-ray binary LS 5039, its possible association 
with SNR G016. 8-01.1, and the motivation to obtain the proper 
motion of the pulsar PSR J1825-1446. In Sect. [3] we describe 
the VLBA astrometric project on PSR J 1825 -1446, and we ob- 
tain the proper motion of the source. In Sect. |4]we present three 
sets of radio observations of LS 5039, as well as the available 
radio and optical astrometry of the source, and we compute the 
proper motion and discuss the uncertainties of the fit. In Sect. [5] 
we compute the Galactic velocity of both sources. In Sect. [6] we 
analyse the past trajectories of the sources and discuss their pos- 
sible birth location. Finally, we summarise the obtained results 
in Sect. 13 



2. Compact objects in the field of SNR G01 6.8-01.1 

2.1. The gamma-ray binary LS 5039 

LS 5039 is a gamma-ray binary system that displays non- 
thermal persistent and variable emission from radio frequen- 
cies to high-energy (HE; E > 100 MeV) and very-high-energy 
(VHE; E > 100 GeV) g amma rays. The s ystem contains 
a brig ht ON6.5 V((f)) star dClark et al.l 120011 McSwain et all 
2004) and a compact object of unknown nature. This degen- 
erate companion has a mass greater than ~ 1.5 M , which 
depends on the b inary system inclin ation, which is an un- 
known parameter (Cas ares et al.l 12.005b . The orbit of the sys- 
tem has a period of 3.9 days and an eccentricity around 0.35 
dCasares et al.1 l2005UAragona et al.ll2009HSartv et al.l20lTl) . The 
distance to the source has recently been updated to 2.9 + 0.8 kpc 
dCasares et al.ll2012al). The persistent sy nchrotron radio emis- 
sion dMartf et all 1 19981; IRibo et al.l [19991; iGodambe etail 120081: 
Bhattacharvva et al.l 120121) appears extended when observed 
Baseline In terferometry (VLBI), at scales of 5-300 m illiarcsec- 



ond (mas) dParedes et al.ll2000l l2002t IRibo et al.ll2008l). The X- 
ray sp ectrum of the system dReig et al.ll2003t iMaitocchi a et alj 
2005) is well fitted by a simple powe r law, and it is clearly peri- 
odic, as shown by RX T E and Suza ku (Bosc h-Ramon et ail 2005: 
Taka hashiet al.ll2009l) . iDurant et all (1201 lh discovered a large- 
scale extended component in X-rays up to 1' from LS 5039. The 



system has also been detected by INTEGR AL up to 200 keV 
dGoldoni et all 12001 iHoffmann et al.l l2009l) . and with BATSE 
up to 1 MeV dHarmon et al.ll2004l). LS 5039 wa s first associated 
with an EGRET source in iParedes et al.l d2000), and variability 
was dete cted by Fermi bet ween 100 and 300 GeV with a 3.9- 
d period (Abd o et al.ll2009l) . The system is also a TeV emitter, 
with persi stent, variable, and pe r iodic emission, as detected by 
H.E.S.S. dAharonian et al.ll2005l 12006). No short-period pulsa- 
tions were fou nd that could demons trate the pre sence of a pulsa r 
either in radio (IMcSwain et al.l201ll) or X-rays dRea et al.l201ll) . 
We note that traditionally the name of the optical star, LS 5039, 
has also been used to refer to the binary system, and so we will 
do in this paper. 

In all gamma-ray binaries, the nature of the compact 
object is fundamental for understanding the physical pro- 
cesses involved in the particle acceleration that is responsi- 
ble for the multi-wavelength emission. If the compact ob- 
ject is a black hole, the accelerated particles would be pow- 
ered by accretion, and produced in the jets o f a microquasar 
dParedes et al.ll2 006; Bosch-Ramon etal||2006j and the review in 
iBosch-Ramon & Khan gulvanl 1 20091) . On the other hand, if the 
compact object is a young non-accreting pulsar, the particle ac- 
celeration would be produced in the shock between the relativis- 
tic wind of the pulsar and the stellar wind of the massive com- 
panion star. This scenario was des cribed in | Maraschi & Treves! 
dl98l|) . iTavani & Aronsl (fl997|) . iKirketalJ d 19991) iDubus 
(2006), Khang ulvan et alJ 00071) . and iBogovalov et al] 12008. 
2012). Dedicated discussions on the wind-wind collis i on sce - 
nario for LS 5039 can be found in iKhangulvan et al] d2007l) . 
Sierpowska-Bartosi k & Torresl (120071) . IDubus et aH (j2008), and 
IKhangulvan et alj d2008l) . Considerations on the spectral sig - 
nature of the pulsar wind are provided in ICerutti et al.l d2 008 ), 



on secondary casc ading in IBosch-Ramon et aT ( 20081) and 



ICerutti et al.l (120101). while X-ray abso rption and occultation are 
discussed in ISzostek & Dubusl d201 ll) . a nd limits on the spi n- 
down luminosity of the putative pulsar in|Zabalza et aT] d201 ll) . 

The orbital parameters of LS 5039 such as the eccentric- 
ity, masses, and orbit i nclination are not well constrained or are 
still under discus sion dCasares et al.ll2005l lAragona et al.l [2009; 
ISartv et al.ll20Tlt) . The proper motion and the sp a ce v elocity of 
the source were first determined in IRibo et all (120021) . These 
authors found that LS 5039 is escaping from its own regional 
standard of rest with a total systemic velocity of 150 km s _1 . 
Th e past trajectory of LS 5039 for the last 10 5 yr computed 
in IRibo et all d2002l) marginally suggests an association with 
SNR G016. 8-01.1. These authors also discovered an Hi cavity 
in the ISM, and argued that it may have been caused by the stel- 
lar wind of LS 5039 or by the progenitor of the compact object in 
the system. An in-depth discussion on the supernova process that 
pr oduced the compact objec t in the system LS 5 39 can be found 
in IMcSwain & Giesl d2002T), IRibo et al.1 d2002l) , IMcSwain et alj 
d2004l) . and lCasares et al.l (120051) . The authors discussed the ex- 
pected mass loss during the supernova explosion that produced 
the compact object, obtained from the orbital parameters and 
the mass of the system. This allowed them to estimate some or- 
bital parameters of the pre-supernova system. We note that the 
main limitation for obtaining accurate pre-supernova properties 
of LS 5039 is the determination of the orbital parameters, which 
can be improved by means of optical observations. 

The region around LS 5039, shown in Fig. [1] contains 
a complex rad i o st ructure that has initially been described 
in IReich et aD dl986l) . By means of observations with the 
Effelsberg 100-m telescope at 1.4, 4.8, and 10 GHz, these 
authors identified the complex as a composition of different 
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Table 1. Parameters of the isolated pulsar PSR J 1825- 1446. 
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Fig.l. Wide field map around SNR G016.8-01.1. The grey 
scale corresponds to the Parkes-MIT-NRAO survey at 5 GHz 
dTasker et al.ll 19941) . and the contours, with higher resolution, to 
the N RAO VLA Sky Survey (NVSS) at 1.4 GHz dCondon et alJ 
1998). The dashed line outlines the area co vered by the diffuse 
emission at 4.75 GHz in lReich et'a l. ( 1986), and the dotted line 
the polarised emission in the same image. These two compo- 
nents trace the morphology of SNR G016. 8-01.1. The more rel- 
evant sources in the field are indicated. 



objects. The brightest source corresponds to the H n region 
RCW 164, which co i ncide s with the optical 8' x 6' counter- 
part in Rodger s~et al.l (Il960l) . It appears to be an optically thick 
thermal component that produces the inverted integrated radio 
emission of the global structure. Other sources within the struc- 
ture are a source with 0.11 Jy at 4.8 G Hz coincident wit h the 
Wolf-Rayet star WR 115 (IC14-19 in lReich et all 1 19861) . and 
an unidentified source with 0.13 Jy at 4. 8 GHz, also detecte d 
with the VLA as NVSS J182535-145555 dCondon et al.lll998h . 
Finally, diffuse emission with a size of ~ 30' surrounds all the 
sources. The source is also visible at 2.7 GHz in the Effelsberg 
11 cm radio continuum survey of the Galactic plane (IReich et al.l 
1 1 990t) . This emission was identified as a supernova remnant be- 
cause it is highly polarised (up to 40% at 4.8 GHz. IReich et alJ 
119861) . The dashed and dotted contours in Fig. [TJ delimit the dif- 
fuse and polarised emission, respectively, and should trace the 
shape of the supernova. Unfortunately, the high contamination 
of the bright RCW 164 prevents the determination of the SNR 
parameters, and consequently its age, distance, centre, and sur- 
face brightness are unknown. The lack of polarization in the di- 
rection of RCW 164 indicates that the SNR is located behin d 
RCW 164, which is located at ~ 1.8 kpc dRibo et alJ 12002). 
Therefore ~ 2 kpc is the lower limit for the SNR distance. 
The size of the polarised structure is ~ 30' and it is approx- 
imately centred at ffnooo o = 1 8 h 25™3 and feooo.o = -14°46' 
dReichetalJll986t [Green 20091) . 



2.2. The nearby isolated pulsar PSR J1825- 1446 

As part of our project to unveil the origin of the binary LS 5039, 
we searched for other possible sources, such as fossil neutron 



r arameter 


oymbol 


Value 


i\.ei. 


Pulsar period 


P 


0.279186875177(5) s 


(1) 


Frequency 


./o 


3.58183026823(6) 


(1) 


Freq. derivative 


fi 


-2.909526(5) x lO" 13 s" 1 


(1) 


Freq. second derivative 


h 


2.24(9) x 10~ 25 s~ 2 


(I) 


Characteristic age 


Tc 


1.95 x 10 s yr 


(1) 


Dispersion measure 


DM 


357 ±5 


(1) 


Rotation measure 


RM 


-899+ 10 rad nr 2 


(1) 


Distance (DM) 


d 


5.0 ± 0.6 kpc 


(2) 


Spindown luminosity 




4.1 x 10 34 ergs" 1 


(1) 


Flux density 


S 0.6GHz 


1.8±0.3mJy 


(3) 




S 1.1GHz 


2.4 ± 0.5 mjy 


(3) 




S 1.4GHz 


3.0 ± 0.2 mjy 


(4) 




S 1.6GHz 


2.5 ± 0.2 mjy 


(4) 




S 4.8GHz 


1.2 ±0.1 mjy 


(3) 


Pulse width (1.4 GHz) 


w w 


27 ms 


(1) 


Pulse width (4.8 GHz) 


w 10 


6.8 ms 


(3) 


Polarisation 




90 ± 6% 


(5) 


Proper motion 


jx a cos 6 


10 ± 18 mas yr~' 


(1) 




H5 


19+115 mas yr" 1 


(1) 



Notes. The proper motion is updated in this paper, see text. 



References. (1) Hobbs et al. 1 200 4]); (2) using t he model in 

ICordes & Laziol d2002l); (3) iKiiak et alJ <2007l) : (4) lLorimer et all 
j!995h ; (5) IVonH oensbroechl dl998l) . 



stars, that may be related to SNR G016. 8-01.1. The aim was 
to obtain the proper motion by means of VLBI of the best can- 
didates to check any possible relation with the SNR. The most 
obvious candidate is PSR J1825-1446, which, in projection, lies 
within the SNR itself. To explore other possi ble candidates we 
used the ATNF Pulsar Catalogue, PSRCAlQ (Manchester et al. 
120051) . which is the most complete and updated pulsar database, 
and accepts searches based on pulsar parameters. We restricted 
the search to those angular distances for which the pulsars proper 
motion and age make their past trajectory compatible with the 
location of the SNR. We note that this selection is not complete 
because, among the known pulsars many do not have good deter- 
minations of the proper motion and/or the age, and their trajec- 
tories are compatible only because of their large uncertainties. 
We discarded pulsars with distances below 1.9 kpc and above 
9.0 kpc. We found no other candidate with enough constraints to 
justify a long-term VLBI project, and therefore we restricted our 
observations to PSR J1825-1446. 

PSR J1825-1446 (B1822-14) is an isolated pulsar with a 
period of 0.28 s discovered in a high radio freque ncy survey 
conducted at Jodrell Bank by Clifton & Lvnd dl986l) . The main 
properties of the pulsar are described in Table [TJ The rad i o spec - 
trum of PSR J 1825 -1446 can be found in Kii ak et alJ d2007l) . 
and it presents a turnover at 1.4 GHz. The spectral index is 
-1.4 ± 0.7 above this frequency, and it is inverted at lower fre- 
quencies. The proper motion in Table[TJwas obt ained by means 
of pul sar timing. The pulse profile obtained by iGould & Lvnel 
shows a single component at 1.4 and 1.6 GHz. We com- 



puted the distance to PSR J1825-1446 using the pulsar disper- 
sion measure, although it is not a direct measurement because it 
depends on the model of the electron distribution in the Galaxy. 
We used the dispersion measure quoted in Table [TJ and the 
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electro n density model NE 2001, described in ICordes & Laziol 
(120021) . and obtained a distance of 5.0 + 0.6 kpc. The distance 
uncertainty is 10%, although it could be overestimated if the 
medium along the line of sight of the pulsar is denser than given 
by the model. 

The possible association between SNR G016. 8-01.1 
and PSR J1825-144 6 has been discussed in the literature. 
IClifton & Lvnel d!986l) already suggested a possible association 
based only on a positional coincidence. Clift on et alj d 19921) 
mentioned the possible relation between them, and the proximity 
with the intense H n region RCW 164 (see Fig. [I}- However, they 
already found that this H n region appears to be much closer to 
the Sun than the pulsar, thus contributing to the DM of the pulsar. 
The authors also argued that the probability of having a pulsar in 
a randomly selected area of the size of SNR G016. 8-01.1 at this 
Galactic longitude is about 1 in 18, and therefore the positional 
proximity may only be by chance. Some notes on associations 
between old (t > 10 s yr) pulsars and supernova remnants can be 
found in lGaensler & Johnston! (Il995l) . 

3. Astrometry and proper motion of 
PSR J1 825-1 446 

The previously known proper motion of the pulsar, shown in 
Table Q] has large uncertainties and did not allow us to study the 
origin of this pulsar. The astrometry in the literature, detailed in 
Sect. I3.2I is scarce. Therefore, we conducted a dedicated 2-yr 
VLBI astrometric project to measure the proper motion of the 
pulsar. 

3.1. Radio VLBI observations 

Before starting the VLBI project, we conducted a short VLA 
observation at 8.6 GHz to obtain the initial position required 
for the correlation of the VLBA data. The VLA observation 
was conducted on November 14, 2008, in the A configuration 
of the VLA, and lasted 30 minutes. At that epoch the VLA 
was performing the upgrade to the EVLA. Seventeen antennas 
used the new EVLA feeds, eight used the old VLA feeds, and 
four antennas were not available. We observed the source dur- 
ing three scans of 3.6 minutes, bracketed by 1.6-minute scans 
on the phase calibrator J1832-1035 (1832-105). The correla- 
tion position of the reference source J1832-1035 was aj2oooo = 
18 h 32 m 20 s :8436and5j2ooo.o = -10°35'11"299. Standard datare- 
duction in AIPS0 was performed, including standard normaliza- 
tion of the flux density between EVLA and non-EVLA antennas. 
The data were imaged with an intermediate robust parameter of 
within AIPS, as a compromise between angular resolution and 
sensitivity, providing a synthesized beam of 460 x 210 mas at 
a position angle (PA.) of 37°. We used a cell size of 10 mas. 
The image shows a point-like source with a peak flux density of 
0.76 + 0.07 mJy beam -1 , and a total flux density of 0.8 + 0.1 mJy. 
The pulsar position was measured with the AIPS task JMFIT. 

The VLBI campaign consisted of three observations with the 
Very Long Baseline Array (VLBA) at 5 GHz separated by one 
year. The observations were conducted on the same day (May 4) 
in 2009, 2010, and 201 1, during six hours between UTC 07:40 
and 13:40. Observations separated by one year prevent any an- 
nual parallax displacement, which is expected to have a maxi- 
mum amplitude of ~ 0.2 mas for a source at 5 kpc. Additionally, 
identical schedules were chosen to minimise any observational 

2 The NRAO Astronomical Image Processing System. 
http://www.aips.nrao.edu/ 



Table 2. Ionospheric correction applied to the VLBA data on 
PSR J1825-1446. Units are mas. 



Epoch 


Mean offset 


Dispersion 


Best 




a 


6 


a 5 




2009 
2010 
2011 


0.13 
-0.05 
0.03 


0.92 
1.96 

2.73 


0.04 0.17 
0.09 0.36 
0.19 0.35 


CODE 
IGS 
IGS 



Notes. CODE: Center for Orbit Determination in Europe. IGS: the 
International GNSS Service 

difference between epochs (e.g. source elevation, scan structure). 
The observations were conducted at a frequency of 5.0 GHz 
(6 cm wavelength), which was chosen as a compromise between 
angular resolution, spectral index of the pulsar (a ~ —1.4), 
and the quality of the best VLBI calibrator next to the source, 
J1825-1718, which is partially resolved at long baselines due to 
interstellar scattering. The first epoch was recorded with a total 
data rate of 256 Mbps per station distributed in four dual po- 
larisation subbands, each of them with a bandwidth of 8 MHz, 
and recorded using 16 frequency channels, two-bit sampling, 
and 2 s of integration time. The correlation of this first epoch 
was performed in the hardware VLBA correlator at Socorro. 
The second and third epochs, observed and correlated in 2010 
and 2011, respectively, were recorded at the available bit rate 
of 512 MHz, and eight subbands were used, which improved 
the sensitivity of the observation. The 2010 and 2011 co rrela- 
tions were performed with the DiFX (De ller et al.l 12007) soft- 
ware correlator at Socorro. The data from the three observations 
of PSR J1825-1446 were also correlated using pulsar gating 
(see Appendix lAl. 

The observations were performed using phase-referencing 
on the nearby phase calibrator J1825-1718, located at 2.5° 
from the pulsar, switching between them with a cycling time 
of four minutes, which is compatible with the expected coher- 
ence time. The phase calibrator was correlated at ffnooo.o = 
18 h 25 m 36 s .53228 and feooo.o = -17°18'4978485. However, for 
the final fit result we used the last available position of the source 
in the frame of ICRF (see Sect. 13 . 3b . As an astrometric check 
source we observed the quasar J 1844- 1324, which is separated 
5° from the pulsar and 6° from the phase calibrator. The source 
has a flux density of 65 mJy, which, in natural weighted images, 
is dominated by a compact core that accounts for 94% of the to- 
tal emission, and displays slightly extended emission eastwards. 
Only the compact core is seen in the uniformly weighted im- 
ages used to measure its position. Two scans of the fringe finder 
J1733-1304 were observed at each epoch. 

The data reduction was performed in AIPS. Flagging based 
on predicted off-source times, owing to slewing or failures, 
was applied using UVFLG. A priori visibility amplitude cal- 
ibration used the antenna gains and the system temperatures 
measured at each station. We used ionospheric total electron 
content (TEC) models based on GPS data obtained from the 
CDDIS data archival to correct the data from the sparse iono- 
spheric variations. Several ionospheric models are produced 
each day by different groups (i.e. the Jet Propulsion Laboratory 
(JPL), the Center for Orbit Determination in Europe (CODE), 
the Geodetic Survey Division of Natural resources Canada 

3 The Crustal Dynamics Data Information System 
http://cddis.nasa.gov/ 
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Table 3. Parameters of the observations of PSR J 1825 -1446, 
and the resulting images. The last column is the improvement 
factor of the signal-to-noise ratio achieved when using pulsar 
gating. 



MJD 


Array 


Freq. 


HPBW 


(S/AOo 


(S/N)o 


Gain 






[GHz] 


[mas 2 at °] 








55607.83 


VLA 


8.6 


460x210 at 37 


11 






54955.44 


VLBA 


5.0 


4.9x2.2 at 5 


3.5 


9.8 


2.8 


55320.44 


VLBA 


5.0 


6.2x2.7 at 6 


6.6 


11.5 


1.7 


55685.44 


VLBA 


5.0 


6.9x2.9 at 1 


6.1 


8.4 


1.4 



(EMR), the ESOC Ionosphere Monitoring Facility (ESA), and 
the Universitat Politecnica de Catalunya (UPC), among others). 
We applied the corrections from all available models to the data 
to check the consistency between models. There are currently 12 
models available for the data obtained in 2009, 16 for 2010, and 
19 for 201 1. In Table [2] we show the offset position between the 
uncorrected data and the corrected data (average value consider- 
ing all models) for right ascension and declination, respectively. 
We also show the standard deviation of the positions measured 
with the different models, and the model that provides a posi- 
tion closer to the mean value, which was used for the final data 
reduction. 

Standard instrumental corrections were applied (parallactic 
angle, instrumental offsets, and bandpass corrections). Fringe- 
fitting on the phase calibrator was performed with the AIPS 
task FRING. The amplitude and phase calibration, flags, and 
bandpass correction tables were applied to the target pulsar and 
the astrometric check source data, which were averaged in fre- 
quency, and clean images were produced with IMAGR. A cell 
size of 0.1 mas was used for cleaning all images, which were 
produced using a weighting scheme with robust parameter -2 
within AIPS. A tapering of 60 MA was applied to avoid the use 
of phases where the correlated flux density of the calibrator was 
below ~40-10 mJy. No self-calibration of the data was possible 
because of the low flux density of the target source. 

The source has a declination of -15°, and the maximum el- 
evation for most of the VLBA antennas is 40°. Some phase de- 
coherence was present when using the data from the beginning 
and the end of the observations, presumably as a consequence 
of the atmospheric effects caused by the low elevation of the 
source and the phase calibrator. To reduce the astrometric er- 
rors, we flagged all visibilities with elevations below 20°, and 
only used the data closer to the culmination of the source, dis- 
carding about 1.5 hours at the beginning and at the end of the 
observations. To reduce the effects of the extended emission of 
the astrometric check source, visibilities with Mv-distances be- 
low 10 MA were not used. The positions of PSR J1825-1446 
and the astrometric check source were measured by fitting a 
Gaussian component with JMFIT within AIPS. In Table [3] we 
show some parameters of the observation, the resolution, and a 
comparison of the obtained signal-to-noise ratio with and with- 
out pulsar gating (see Appendix lAl for the details). The position 
errors of PSR J 1825 -1446 were finally computed as the stan- 
dard deviation of the three positions of the astrometric check 
source. The final positions measured in the images from the data 
correlated with pulsar gating are shown in Table |4] 



3.2. Archival astrometry of PSR J1825- 1446 

To measure the proper motion of PSR J 1825 -1446 we 
sea rched all currently available ast rometry of this isolated pul- 
sar. |^rafl&_Sc^OTirighause^(lT997|) imaged the field around sev- 
eral pulsars with the VLA. They observed PSR J 1825 -1446 
(B1822-14) on March 14, 1993 at 8.4 GHz for two scans 
of 10 minutes each. The observations were conducted with 
the VLA in its B configuration. The phase calibrator was 
J1832-1035 (J1832-105, 1829-106). We have re-analysed the 
data to measure the position of the pulsar. Although the scan 
duration is longer than optimal for astrometry, the pulsar is 
detected as a point-like source with a peak flux density of 
0.78 ± 0.10 mJy beam -1 when using uniform weighting, which 
yields a synthesized beam size of 1.0 x 0.6 arcsec 2 . We trans- 
formed the positions of the observation from B1950 to J2000, 
and we updated the position of the phase calibrator to the best 
currently known position, quoted in Sect. 13.11 We fitted the po- 
sition using JMFIT within AIPS, and the errors correspond to 
formal errors of a Gaussian fit. The fitted position can be found 

in Table |4] 

The position of PSR J 1825 -1446 found in iHobbs et all 

(2004), obtained by means of pulsar timing, is also shown in 
Table [4] The position is in the solar system centre of mass, or 
barycentric, frame, and it was obtained by fitting the pulsar pe- 
riod and period derivative, the ecliptic position and proper mo- 
tion, among other variables, to the time of arrival (TOA) of the 
pulses during several years. However, we can see in Table 2] 
that the timing position differs from the expected position by 
245 mas (3-cr) in right ascension, which appears to indicate that 
the uncertainty may be underestimated (see also Fig. [2] and the 
difference with the declination error). The difference between 
the errors in right ascension and declination can be understood 
considering the pulsar position. The pulsar is located near the 
ecliptic plane (with ecliptic latitude ft = 8°53), which has two 
consequences in the errors of the ecliptic coordinates (A, ft). In 
particular, the uncertainty in A and ft are inversely proportional 
to cos ft, and sin/3, respectiveljQ. For PSR J1825-1446 we have 
that (cos/Jr 1 ~ 1.0, and (sin/JT 1 ~ 6.7, which intrinsically 
yields larger errors in the determination of the ecliptic latitude 
by means of pulsar timing. Large errors in the ecliptic latitude ft 
produce the large error in declination, which is indeed a factor 
6.7 larger than the uncertainty in right ascension. A nother possi- 
ble sou rce of error is the pulsar timing fit solution. IHobbs et al.l 
(120041) fitted 197 time of arrival (TOA) data points spanning 
17.4 years and obtained a proper motion of fi a cos6 = 10+18 
and fig = 19 ± 115 mas yr _1 . The fitting time residuals were 
7.539 ms. If this value were understimated, this would explain 
the small error in right ascension. Another possible source of 
error is the uncertainty in the link between the barycentric celes- 
tial reference system (BCRS) and the ICRF. However, it should 
be well below the 245 mas difference found between the tim- 
ing position_and_fhe exp ected position at that epoch; for exam- 
ple IChatteriee et al.l (120091) found differences between both ref- 
erence frames below 1 mas. We have not found any clear ex- 
planation to justify the discrepancy between the right ascension 
position given by pulsar timing and the expected one. Given the 
large difference, and the intrinsic difference between an interfer- 
ometric measurement and a timing fit, we have not included the 
timing position in our final determination of the proper motion 
of PSR J1825-1446, although we have included the position in 
Tableland Fig. [2] for reference. 



http://www.cv.nrao.edu/course/astr534/PulsarTiming.html 
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Table 4. Astrometry and fit residuals of PSR J 1825 -1446. 



MJD 


Epoch 


Array 


Freq. 


Q'nooo 


CT qi cos 6 


<5j2000 




R.CS.Q. cos s 


Res. ,s 




Y-M-D 




[GHz] 


[h, m, s] 


[mas] 


[°, "] 


[mas] 


[mas] 


[mas] 


49060.53 


1993-03-14 


VLA-B 


8.4 


18 25 2.94(3) 


40 


-14 46 52.77(6) 


60 


-37 


-61 


49480.00 


1994-05-08 


Timing 


0.4-1.6 


18 25 2.927(5) 


75 


-14 46 52.6(5) 


500 


-245 


142 


54784.04 


2008-11-14 


VLA-A 


8.6 


18 25 2.9539(8) 


12 


-14 46 53.163(15) 


15 


13.4 


-0.0 


54955.44 


2009-05-04 


VLBA 


5.0 


18 25 2.95324(4) 


0.6° 


-14 46 53.1767(5) 


0.5* 


-0.6 


-0.3 


55320.44 


2010-05-04 


VLBA 


5.0 


18 25 2.95398(4) 


0.6° 


-14 46 53.2047(5) 


0.5* 


0.7 


0.7 


55685.44 


2011-05-04 


VLBA 


5.0 


18 25 2.95460(4) 


0.6° 


-14 46 53.2347(5) 


0.5* 


0.0 


-0.4 



Notes. (fl) An additional uncertainty of 0.9 mas in right ascension due to the global astrometry errors on the phase-reference source position should 
be included when considering this position in the ICRF. '*' The phase-reference source has an additional uncertainty of 1.4 mas in declination in 
the ICRF. 
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Fig. 2. Position of PSR J1825-1446 with respect to time (MJD in lower axis and year in top axis) in right ascension (left) and 
declination (right). The solid black lines represent the fitted proper motion, and the grey area the uncertainty of the fit at 3-cr level. 
The black squares represent the VLA positions, the black circles indicate the VLBA positions, and the grey circle the position from 
pulsar timing (not used in the fit). The inset shows the position residuals from the fitted model for the 2008-2012 period. The source 
positions are plotted with uncertainties at l-cr level. 



3.3. Proper motion of PSR J 1825- 1446 

To obtain the proper motion, we used the two VLA and the 
three VLBA observations presented in Sect. [3] The reference 
position of the phase calibrator for the VLA observations was 
updated to the currently known position of the phase calibra- 
tor, quoted in Sect. I3.ll All VLBA positions were updated ac- 
cording to the most recent available position^ of J 1825- 17 18 
arjaooo.o = 18 h 25 m 36^532278 ± (£00006 (or ±0.9 mas) and 
<5j2(>oo.o = -17°18'49'.'8482 ± 070014 (or ±1.4 mas). This global 
uncertainty of the reference source position should be included 
when working in the ICRF. However, in this case it is only a 
common offset shared by all VLBA observations. We have com- 
puted the proper motion fit using the position of J 1825 -17 18 
as a reference. This means that the 0.9 and 1 .4 mas uncertain- 
ties have been included in all the non-VLBA (but not in the 
VLBA ones), to compute the proper motion, and it has been 
added afterwards to the final fitted position. All positions and 
uncertainties are shown in Table |4] We fitted the proper mo- 

5 GSFC 2010a astro solution, from sched vlO.O 



tion of the source as a linear model of position as a function 
of time for right ascension and declination separately. We fitted 
the data using a L evenberg-Marquardt nonlinear least-squares fit 
dMarquardtl 19631) . The solution is highly dominated by the three 
VLBA observations, which were performed in the same sched- 
uled conditions (see Sect. 13. II ). The reference epoch is computed 
as the weighted mean MJD, which is MJD 55291.525. The ob- 
tained astrometric parameters for PSR J 1825 -1446 are 

ffj20oo = 18 h 25 m 295389 (±0.3 ± 0.9 mas) , 
<Sj20oo = -14° 46' 53V2031 (±0.3 ± 1.4 mas) , 
H a cos 8- 10.0 ± 0.3 mas yr , 
Us = -29.0 ±0.3 masyr" 1 , 

where the first set of errors in the positions corresponds to statis- 
tical uncertainties, and the second set corresponds to the uncer- 
tainty in the phase reference position. The degrees of freedom of 
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Table 5. Parameters of the LS 5039 observations described in 
Sect. 3] parameters of the resulting images, and the total and 
peak flux densities measured. 



MJD 
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Freq. 
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S. 




S peak 






[GHz] 


[mas 2 at °] 


[mJy] 


[mjytr 1 ] 


52856.20 


VLA+PT 


15.0 


148x39 at 23 


9.6 ± 


0.5 


9.0 ± 0.3 


53167.35 


VLBA 


8.5 


4.5x2.2 at 7 


10.9 ± 


1.5 


6.0 + 0.6 


53532.35 


VLBA 


8.4 


2.8x1.1 at 6 


7.1 ± 


0.9 


2.6 + 0.3 


53764.71 


VLBA 


8.4 


3.2x1.0 at 10 


4.4 ± 


0.7 


3.1 +0.3 


55635.15 


EVN 


8.4 


5.6x1.3 at -15 


8.0 ± 


2.0 


6.1 ± 1.0 



the fit are 4, and the reduced x 2 is 0.95 and 0.98 in right ascen- 
sion and declination, respectively. 

4. Astrometry and proper motion of LS 5039 

To fit the proper motion of LS 5039 we used new interferomet- 
ric radio data together with archival radio and optical data. In 
this section we will only describe in detail the data reduction of 
three radio projects. The description of the other data used is al- 
ready published, or will be published elsewhere. In particular we 
will describe a VLA+PT observation conducted in 2003, three 
VLBA observations conducted between 2004 and 2006, and a 
short EVN observation conducted in 2011. The parameters of 
these observations are summarised in Table [5] 

4. 1 . VLA-PT observation 

We observed LS 5039 with the VLA plus the VLBA antenna 
in Pie Town (PT) at 15.0 GHz on August 05, 2003. We used 
the VLA in its A configuration plus PT, which is located at a 
distance of 52 km from the VLA site, and provides longer base- 
lines with all the VLA antennas. The observation spanned from 
01:00 to 08:40 UTC, and was centred at MJD 52856.2. The 
phase reference calibrator was J1832-1035 (1832-105), located 
at 4.5° from LS 5039. Two IFs of 50 MHz were used. The cy- 
cle time including the phase calibrator and the target source was 
3 minutes. The ICRF source J191 1-2006 (1911-201) was ob- 
served regularly during the observations. The flux calibrator was 
3C 286 (1331+305). The orbital phase of the binary system at 
the time of the obse rvation was 0.92, using the ephemeris from 
ICasares et al. I (120123) . 

Standard reduction procedures were applied to the data, and 
a precise position of the source was obtained when using a uni- 
form weighting scheme (robust -5 within AIPS), and a cell- 
size parameter of 5 mas. The source position, fitted with JMFIT 
within AIPS, has an uncertainty of 1 and 2 mas in right ascen- 
sion and declination, respectively. The observational details, the 
obtained synthesized beam size, and the fitted components are 
shown in Table [5] 

4.2. VLBA astrometric project 

The VLBA projects BD087 and BD105 are part of a long-term 
astrometric project to measure the proper motion and parallax of 
several X-ray binaries. The segment BD087G, and the continua- 
tions BD105A and BD105G observed LS 5039 at a frequency of 
~8.5 GHz. The observations took place on June 11, 2004, June 
11, 2005, and January 29, 2006, and lasted 6, 6, and 5 hours, 



respectively. The three observations were recorded with a to- 
tal data bit rate of 128 Mbps. In the first epoch, the data rate 
was distributed in 8 IF channels with a bandwidth of 4 MHz 
each, with one single polarisation (right circular polarisation) 
and two-bit sampling, while the latter two used two IF chan- 
nels with 8 MHz bandwidth and dual circular polarisation, also 
recorded with two-bit sampling. 

The phase reference source was J 1825 -17 18, which is 
slightly resolved also at 8.4 GHz (although it is much more 
compact than at 5 GHz), and is located at 2.47° from LS 5039. 
The cycling time of the phase-referenced observations was three 
minutes. Several check sources were observed during the three 
epochs. In particular, the sources N1834-1433, N1819-1419, 
and N1818— 1 108 were observed in all epochs and were used 
to estimate the astrometric uncertainties. The fringe finder was 
J1733-1304. The data reduction is analogous to the one de- 
scribed in Sect. 13.11 In this case, an statistical approach to de- 
termine the best ionospheric correction was not applied, because 
the effect is smaller at this frequency. We note that the correction 
of the Earth orientation parameters was not applied to the data 
because it did not affect the position measurement, and there- 
fore we preferred not to apply this additional correction to the 
datasets. 

Owing to the lower flux density of LS 5039 at 8.4 GHz, the 
three epochs were imaged using a weighting scheme with robust 
parameter within AIPS, as a compromise between spatial reso- 
lution and signal-to-noise ratio. To reduce the effect of the source 
extended emission (see Sect. 14.51 ) we omitted the baselines with 
shortest wv-distances. We found a compromise between signifi- 
cance of the detection, and the wv-distance reduction by cutting 
the baselines below 20, 30, and 40 MA for the three epochs. For 
the first epoch a lower threshold was required because the data 
quality was worse. On the other hand, the longest baselines had 
to be softened because of the poor quality of the partially re- 
solved phase reference source. We reduced the weight of base- 
lines longer than a certain value. For the first epoch we used 
a tapering parameter of 60 MA. For the latter two projects, the 
tapering was at 120 MA. The resulting resolution and fitted com- 
ponents are shown in Table [5] 

4.3. EVN observations 

The most recent VLBI observation of LS 5039 was performed 
with the European VLBI Network (EVN). We observed several 
sources around LS 5039 to search for possible phase calibra- 
tors and in-beam calibrators, and test the longest EVN baselines 
to check the viability of an astrometric project. Among other 
sources, we observed the phase calibrator J1825-1718, also 
used in the previously described VLBI observations. The obser- 
vations were conducted on March 15, 201 1, between UTC 02:00 
and 05:00. We only observed 11 scans on LS 5039, brack- 
eted by scans on the phase calibrator with a cycling time of 
3.4 minutes and distributed along the observation. The total 
time on-source was approximately 20 minutes. The antennas in- 
cluded in the observation were Effelsberg, Westerbork, Onsala, 
Medicina, Yebes, Svetloe, Badary, Zelenchukskaya, Nanshan 
(Urumqi), Sheshan (Shanghai), and Hartebeesthoek. The data 
were recorded with a total data bit rate of 1024 Mbps, pro- 
viding a high-sensitivity array. The correlated parameters are 
eight IF channels with both circular polarisations, each of them 
with 16 MHz bandwidth (provided by 32 frequency channels of 
500 kHz), and two-bit sampling. 

The data reduction is analogous to the one described in 
Sect. [3T] and Sect. [4/2] No useful FRING solutions were found 
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between Hartebeesthoek and the rest of the antennas, which con- 
siderably reduced the North-South array extend, and therefore 
the resolution in declination on the final image. After the phase 
calibration from J1825-1718, the North-South maximum base- 
line was of about 20-40 MA, while in the East-West direction it 
was of 160-180M/1, although large gaps were present in the uv- 
plane due to the short observation time. We note that large phase 
instabilities appeared due to the low declination of the source. 
We finally imaged the source with robust parameter 0, minimum 
Mv-distance set to 15 MA, and Mv-taper set to 120 MA, Even with 
the long maximum baselines, the final synthesized beam size is 
5.6x1.3 mas 2 in P.A. of -15°. The resolution of the images and 
the fitted components are shown in Table [5] 

4.4. Archival data of LS 5039 

Astrometry of LS 5039 can be obtained from two types of ob- 
servations. On one hand we have positions from optical obser- 
vations of the main star of the system. This observations, com- 
ing from large (and global) astrometric projects, cover a long 
time, of the order of decades, although they provide limited as- 
trometric accuracy of about 10-300 mas. On the other hand, ra- 
dio interferometric observations provide observations of the ra- 
dio nebula around the binary system. These observations have 
been obtained in the last ~10 years, but their uncertainties are 
between 1 and 10 mas. The time span of the optical observations 
can be combined with the precise radio positions to fit an accu- 
rate position and proper motion of the source. However, com- 
bining datasets with considerable different properties brings im- 
portant caveats that have to be taken into account. These will be 
described in Sect. 14.61 The astrometry described here is summa- 
rized in Table. [6] The corresp onding orbital phases were com- 
puted using the ephemerides in Casare s~et al.1 (1201 2al) . 

4.4.1. Archival optical astrometry of LS 5039 

Optical observations of LS 5039 are available from historical 
photographic catalogues since 1905. Updated global solutions 
have been obtained for the historical photographic plates by us- 
ing more recent reference catalogues. The updated versions of 
the classical catalogues provide improved systematic uncertain- 
ties on the star position. However, some newer versions com- 
bine catalogues from different epochs to fit a global solution that 
considers the proper motion of the stars. Therefore, they are not 
suitable for our proper motion determination because they are 
not independent, and they include a proper motion estimation 
in their own position fitting. This is the case, for example, of 
the Palomar Observator y Sky Survey plates used to produce the 
USNO-A2.0 catalogue (Monet] [T998l) . which uses only the blue 
and red origin al plates, while th e new version of the catalogue, 
USNO-B1.0 dMonet et alJl2003l) . with better uncertainties, uses 
data from 1949 to 2002, combined using a proper motion of 4, 
and -14 mas yr~' in right ascension and declination, respec- 
tively. It is also the case of GSC 1 .2 compared to the newer ver- 
sions (GSC 2.2 and GSC 2.3.2), or UCAC1 when compared to 
UCAC2 and UCAC3, which have improved references and sys- 
tematic error treatment. Therefore, although new catalogues are 
currently available, we have used only the versions that are not 
compilation catalogues. This means that most of the optical val- 
ues used here we re already available at the time of writing of 
iRibo etalJd2002l) . 

In particular, for the proper motion determination we used 
the following catalogues. AC2000.2 is the second version of the 



Astrographic Catalog 2000 and provides accurate optical po- 
sitions of LS 5039 on 1906 in Hipparcos Celestial Reference 
System (HCRS, coincident with J2000.0). USNO-A2.0 contains 
a global reduction of the data in the Palomar Observatory Sky 
Survey plates obtained around 1950. The Twin Astrographic 
Catalog version 2 (TAC 2.0) provides a position based on 
four observations (the plates were c ompiled between 1977 and 
1986). Following IRibo et all d2002f) . we included an additional 
uncertainty of 15 ma s in each coordinate, as suggested in 
IZacharias et all (Q~999). We also used the Guide Star Catalog 
(GSC) version 1.2, which consists of a single-epoch collection 
of 6.4° x 6.4° Schmidt plates (the newer versions of the catalog 
are compiled, and proper motions are fitted using Tycho posi- 
tions). The GSC position is less relevant for the fit because of 
the relatively large errors (300 mas) for a position obtained in 
1986. The Tycho-2 Catalog was obtained from the Tycho star 
mapp er observations of the ESA Hipparcos satellite (H0 g et al.l 
2000). It provides the mean position at epoch J2000, corrected 
from proper motion, as well as the obse rved position, wh ich is 
the one used in our analysis. Following Ribo et al. (2002), we 
have used the conser vative error estimate of cr a cos 5 — 149 mas 
and o-g = 120 mas. ISkrutskie etafl (120061) detected LS 5039 
at near-infrared bands with the all-sky Two Micron All-Sky 
Survey (2MASS), and a very precise position was obtained. The 
three observed bands (J, H, and K s ) were tied with the Tycho- 
2 Reference Catalog, so it is in the ICRS frame. Finally, ac- 
curate astrometry of LS 5039 can be found in the latest ver- 
sions of the US Nava l Observatory CCD Astrograph Catalog 
(Zacharias et al H owever, we can only use the position 

from the UCAC1 version (IZacharias et al.l l2000). All these opti- 
cal positions and their uncertainties are quoted in Table |6] 

4.4.2. Archival radio astrometry of LS 5039 

The binary sys t em L S 5039 is a known radio emitter discovered 
bv lMartfetal1(ll998h . Since then, several radio interferometric 
observations have been conducted to study its morphological and 
spectral properties at different orbital phases. Here we list the 
radio projects in which the system has been observed with suf- 
ficient resolution to provide an accurate position of the source. 
In particular, it has been observed with the VLA, the VLBA and 
the EVN at frequencies between 5 and 15 GHz. Observations at 
lower frequencies with the VLA are excluded because they pro- 
vide positions with uncertainties above 07 1 for relatively recent 
epochs, and do not contribute to the proper motion fit when com- 
pared to the 1-10 mas accuracy of the other radio observations. 

We briefly su mmarise the observations included in the fit. 
Mart f et al.l(l 19981) conducted a multi-epoch and multi-frequency 
project of the system with the VLA at four frequencies between 
1 and 15 GHz. They provided a single position derived from 
the 15 GHz data obtained with the highest resolution configura- 
tion (A) of the VLA. The reference source was the ICRF source 
J191 1-2006 (191 1-201), and their error estimate is 10 mas. The 
source was next observed with t he VLBA and the VLA for two 
epochs in June 2000 at 5 GHz dRibo et al]l2008l) . Only one of 
the two epochs provides good astrometry. The source was also 
observed with the EVN within project EF018 in March 2007. 
The project consists of three observations at 5 GHz, observed 
every two days. The low elevation of LS 5039 during the ob- 
servations affects the phase calibration of the data and there- 
fore the astrometry. To obtain a reliable position, we removed 
the data with the lowest elevations and then averaged the three 
measurements. The uncertainty of this position was computed as 
the standard deviation of the three runs. We conducted a similar 
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Table 6. Astrometry and fit residuals of LS 5039 for optical and radio projects. 



Project 


\ JIT** 


Epoch 


0orb 


Q"J2000 




@~a cos <5 


<5j2000 




0~6 


RCS .q. cos 6 


Res. ,j 






Y-M-D 




[h, 


m, : 


>J 


[mas] 


[°, 


"] 




[mas] 


[mas] 


[mas] 


AC2000.2 (O) 


17408.47 


1906-07-17 


0.8(2) 


18 


26 


15.015(8) 


120 


-14 


50 


53.27(11) 


110 


40 


160 


USNO-A2.0 (O) 


33857.61 


1951-07-30 


0.0(1) 


18 


26 


15.034(17) 


260 


— 14 


50 


53.59(25) 


250 


10 


230 


TAC2.0 (O) 


44050.66 


1979-06-26 


0.57(5) 


lb 


Zo 


15.056(5) 


80 


1 A 

— 14 


rn 


54.08(5) 


50 


140 


-8 


GSC1.2 (O) 


46669.35 


1986-08-27 


0.99(4) 


1 Q 
1 K 


ZD 


15.06(2) 


300 


1 A 
— 14 




54.3(3) 


300 


80 


-145 


Tycho-2 (O) 


48559.95 


1991-10-30 


0.01(3) 


1 o 


ZO 


15.043(10) 


150 


1 A 
— 14 


j\) 


54.23(12) 


120 


-150 


-55 


AP357 (R) 


50913.00 


1998-04-10 


0.43(2) 


1 o 


7ft 


15.0559(7) 


10 


— 14 

14 


SO 


54.240(10) 


10 


5.5 


-8.1 


2MASS (0) 


51401.02 


1999-08-11 


0.37(2) 


1 8 

1 o 


26 


15.056(4) 


60 


— 14 




54.23(6) 


60 


-7.3 


18 


UCAC1 (0) 


51649.00 


2000-04-15 


0.85(2) 


18 


26 


15.0563(9) 


13 


-14 


50 


54.277(13) 


13 


-2.9 


-28 


GR021A (R) 


51698.40 


2000-06-03 


0.49(2) 


18 


26 


15.05645(10) 


1.5" 


-14 


50 


54.2512(15) 


2 i, 


-1.6 


-0.5 


Ar4oJ (K) 


jZajb.ZU 


zUUi-Ua-lu 


V.yZ(Z) 


18 


26 


ij.lDBiz(l i ) 


1 .0 


-14 


50 


j4.2sU(2) 


_ 


1 .U 


— 1.0 


BD087G (R) 


53167.35 


2004-06-11 


0.58(2) 


18 


26 


15.05831(14) 


2" 


-14 


50 


54.2865(20) 


2* 


-2.3 


-0.7 


BD105A (R) 


53532.35 


2005-06-11 


0.02(2) 


18 


26 


15.05890(8) 


1.3" 


-14 


50 


54.2964(14) 


1.4* 


-0.4 


-1.7 


BD105G (R) 


53764.71 


2006-01-29 


0.51(2) 


18 


26 


15.05916(7) 


1.1" 


-14 


50 


54.2996(12) 


1.2* 


-1.1 


0.6 


EF018A-C (R) 


54160.71 


2007-03-01 


0.89(2) 


18 


26 


15.0597(2) 


4" 


-14 


50 


54.3092(12) 


1.2* 


-0.2 


0.5 


BR127A-E (R) 


54290.27 


2007-07-09 


0.06(2) 


18 


26 


15.06003(8) 


1.3" 


-14 


50 


54.3096(12) 


1.2* 


1.8 


3.3 


EM085 (R) 


55635.15 


2011-03-15 


0.37(2) 


18 


26 


15.06184(8) 


1.3" 


-14 


50 


54.3452(20) 


2.0* 


2.7 


-0.2 


BC196 (R) 


55803.13 


2011-08-30 


0.37(2) 


18 


26 


15.06188(7) 


LP 


-14 


50 


54.3511(14) 


1.4* 


1.4 


-2.0 



Notes. (o) An additional uncertainty of 0.9 mas in right ascension due to the global astrometry errors on the phase-reference source position should 
be included when considering this position in the ICRF. '*' The phase-reference source has an additional uncertainty of 1.4 mas in declination in 
the ICRF. 



project, observing the source with the VLBA du ring five con- 
secuti ve days in July 2007 at 5 GHz (see details in Moldon et al. 
|2012|) . A similar procedure was used in this case, computing the 
average position of the source from the five images produced 
with uniform weighting. Finally, we used a recent position ob- 
tained within the VLBA large project BC19rJ] (PI: J. Condon), 
whose solutions are av ailable online. The p ositions, uncertain- 
ties (see the overview in lCondon et al.ll201 ll) . and residuals with 
respect to the final proper motion fit are shown in Table |6] 

4.5. Uncertainties in the position of LS 5039 

LS 5039 is known to be extended at mas scales, and a bipolar 
jet-like structure is present in all VLBI (VLBA and EVN) radio 
observations of the source dParedes et al.ll2000l I2002L [Ribo et al. 
2008; Moldon et al] |201 ldl2012ir Although the source displays 
a bright core, we do not know the offset between the peak of 
the radio core and the position of the binary system. There are 
two main reasons. First, the particles accelerated in the system 
travel distances much longer than the system orbit, and the lo- 
cation of the peak of the emission at any moment depends on 
the physical conditions (basically the electron distribution, and 
the synchrotron cooling time). Second, for each observed fre- 
quency, each interferometer array is sensitive to different ranges 
of spatial scales, therefore inter-comparison is not straightfor- 
ward because we are seeing different regions of the particle flow. 
In summary, the unknown position of the system within the ra- 
dio outflow is variable (at time scales of 1 or more days), and de- 
pends on the observed frequency and array. We tried to minimise 
this effect by measuring positions derived from the uniformly 
weighted images, which trace the most compact emission from 
the core component of the source. 

Preliminary fits of the position and proper motion of LS 5039 
using the formal astrometric errors of the radio positions yielded 
reduced^- 2 of 4.1 and 2.3 in right ascension and declination, 

6 https://science.nrao.edu/observing/largeproposals/blackholesearch/ 



respectively, which indicated that a linear model did not com- 
pletely describe the source positions. This means that the for- 
mal astrometric errors for the radio observations underesti- 
mate our real knowledge of the position of the binary system 
within the radio emission. A shift on the peak of the emis- 
sion was observe d in LS I +61 303 PSR B 1259-63, a nd 
HESS J0632+057 dDhawan et alJl2006t iMoldon et al.ll201 1 alibi) 
Based on these observations, and the estimat ions in |bubus 
( 20061) and lBosch-Ramon & Khangulvanl ( 1201 ll) . we expect the 
peak of the emission to follow an ellipse with semimajor axis of 
about ~1 mas, although the size should depend on the frequency, 
and should be larger at low frequencies. It is not possible to mea- 
sure this unknown deviation with the current resolution and data 
sampling, especially when observing at different frequencies and 
resolutions. Therefore, we increased the position uncertainty of 
the radio positions to compensate our ignorance of this time- 
dependent offset. In particular, we added in quadrature an ad-hoc 
uncertainty of 1.0 mas to all radio positions before performing 
the fit. This uncertainty is already included in Table [6] 

Another source of uncertainty is the systematic errors on 
the position of the phase calibrators, which are the reference 
of the measured positions. Fortunately, the phase reference cal- 
ibrator for all VLBI observations was J1825-1718. Although 
the calibrator was correlated at different positions, we have 
shifted all measurements to the common reference position 
quoted in Sect. 13.31 We followed the same procedure as with 
PSR J1825-1446, and these uncertainties were not included in 
the proper motion determination, but were added afterwards to 
the final fitted position. 

We also included the effect of the annual parallax of 
LS 5039, assuming a distance of 2.9 kpc, which corresponds 
to a deviation of 0.3 mas. The equations govern ing the paral- 
lax di splac ements can be found , for example, in B risken et al.l 
(2002J) and lLoinard et al.l (120071) . The additional displacements 
were added to the data before performing the fit, and are not in- 
cluded in Table [6] 
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4.6. Proper motion of LS 5039 

We fitted the reference position and the proper motion to the 
data in right ascension and declination independently for the 
16 positions available, using the Levenberg-Marquardt nonlinear 
least-squares fit. The reference epoch for the fit, computed as the 
weighted mean of the dates of observations, is MJD 53797.064. 
The resulting astrometric parameters for LS 5039 are 

aj20oo = 18 h 26 m 15 s .05927 (±0.5 + 0.9 mas) , 
<5j2000 = -14° 50' 54'.'3010 (±0.5 ± 1.4 mas) , 
fi a cos (5 = 7.10 ± 0.13 mas yr~ , 
fig = -8.75 ± 0.16 mas yr -1 , 

where the first set of errors in the positions corresponds to sta- 
tistical uncertainties, and the second set corresponds to the un- 
certainty in the phase reference position. The reduced x 2 of the 
right ascension and declination fits were 0.75 and 1.35, respec- 
tively. The positions, fit, and uncertainties are shown in Fig. [3] 
The solution is dominated by the radio observations, but the op- 
tical observations provide better stability to the fit. We note that 
the effect of the parallax correction on the data is small. If the 
correction is not included, the reduced \ 2 m right ascension in- 
creases by 0.02, and does not change in declination. We also 
tried to fit the parallax to the data, with and without the ad hoc 
1 mas uncertainty, but the parallax is too small to be measured 
with the current data. 

To better understand the systematic fit uncertainties, we anal- 
ysed the distribution of fi t solutions following a p rocedure simi- 
lar to the one discussed in Chatteri ee et al.l d2009). The bootstrap 
method used consists of fitting the proper motion and reference 
position from different samples of the available positions. For 
each iteration we randomly selected 16 positions, with replace- 
ment, and produced a proper motion fit. By repeating the fit for 
different combinations of positions, we obtained a probability 
distribution of proper motions obtained with all possible combi- 
nations of the data. The probability distribution functions (PDF) 
and the cumulative distribution functions (CDF) of these PDFs 
are plotted in Fig. |4] The distributions were computed with a 
bin width of 0.01 mas yr~'. The number of iterations (different 
combinations with replacement) is 2 x 10 7 , which is more than 
necessary to smoothly cover the one-dimensional range of val- 
ues, but it is required for the 2D distribution (see below). The 
shaded grey areas in Fig. [4] indicate the minimum regions con- 
taining the 68.2, 95.4 and 99.7% of the samples, from darker to 
lighter, respectively. The best-fit value of the proper motion can 
be computed in two ways: as the mode of the distribution (value 
that occurs most frequently), which is marked as the vertical line 
in the top panels, or as the value that accumulates the 50% of the 
events (the median), which is marked as the vertical line in the 
bottom panels of Fig. [4] These two values are not identical, but 
the difference is well below the 1-cr uncertainty of the proper 
motion estimation. We can clearly see that the distribution is ap- 
proximately Gaussian at the central region, but it shows consid- 
erable asymmetric wings. The mode and the confidence intervals 
are fi a cosd - 7.09*2' [j* mas yr~' and fig = -8.82^]' 27 2 mas yr -1 . 

Owing to the asymmetry in the distributions of Fig. [4] we 
computed the 2D distribution of proper motions obtained with 
the bootstrap method in right ascension and declination simul- 
taneously to better identify any other correlation or asymmetry 
between the confidence regions. The simulated fit results were 
gridded with a bin size of 0.01 mas yr~' in both directions. 



The minimum regions covering the 68.2, 95.4 and 99.7% of the 
events were obtained, and they were smoothed with a Gaussian 
filter with a size of 18 bins. The corresponding confidence re- 
gions are shown in Fig. [5] We can see that the distribution is 
clearly skewed towards lower values of fig and more moderately 
to lower values of fi a cos 5. 



8 7 6 -9-8-7 




fi a cos(cS) [mas yr 1 ] n$ [mas yr 1 ] 

Fig. 4. Top: distribution of proper motion fit results for LS 5039 
estimated with the bootstrap method, which is equivalent to an 
un-normalised probability distribution function (PDF) for right 
ascension and declination independently. Bottom: corresponding 
cumulative distribution functions (CDF). For all panels, the grey 
shaded areas correspond to the minimum region that includes 
the 68.2, 95.4 and 99.7% of the events, from dark grey to light 
grey. The vertical line marks the most common value (the mode) 
for the PDF distributions, and the value covering the 50% of the 
events for the CDF. 



5. Galactic space velocity 

We can now study the Galactic space vel ocity of LS 503 9 
and PSR J 1825 -144 6, fol lowing the work in [Ribo et al.l d2002l) , 
Dha wan et alj d2007l) . and lMiller- Jones et al.l (l2009al) . The space 
velocity of each source can be obtained by combining the proper 
motion and the radial velocity of the source (optical systemic 
velocity, y, in velocity lightcurves) to form a space velocity vec- 
tor (fi a cos 6, fig, y). This velocity can be rotated to the carte- 
sian coordinates U, V, W, defined as U positive towards the 
Galactic Centre, V positive towards I = 90° (sense of the 
Galactic rotation), and W positive towa rds the North Galactic 
Pole (N GC). The rotation is described in iJohnson & Soderbloml 
(1987), although we used t he updated definition of th e North 
Galactic Pole direction from Reid & Brunthalerl (12004 (the di- 
rection of the NGP is a 1 ™ 00 = 12 h 51 m 26 s .282 and 4 2000 = 
+27°07'42'.'01, with the zero longitude at 0j 2000 = 122°932). 
We computed the heliocentric Galactic space velocity by sub- 
tracting the peculiar velocity of the Sun with respect to the lo- 
cal s tandard of rest (LSR). We used the recent determination 
from lSchonrich et al.l(l2010l) of U G = ll.lt°f 5 , V = 12.24+° ^, 
W e = 7.25+°'^ km s _1 . We note that V is significantly larger 
than in previous estimations. The obtained heliocentric Galactic 
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Fig. 3. Offsets in right ascension (left) and declination (right) of LS 5039 with respect to time in MJD and year. The solid black lines 
in each panel represent the fitted proper motion, and the grey area the uncertainty of the fit at 3- cr level. The upper panels include all 
data points used for the fit, and the lower panels are the zoom of the last 16 years (marked by the small rectangles), and include the 
residuals with respect to the fitted model. The empty blue triangles represent the optical data, the green squares show the VLA data, 
the red triangles correspond to the VLBI observations at 5 GHz, and the black circles are the VLBI data at 8.6 GHz. The source 
positions are plotted with uncertainties at 1-cr level. 



velocities (U,V,W) LS r for LS 5039 and PSR J1825-1446 are 
shown in Table. [7] 

We note that the distance to the source is the highest 
source of error when converting from angular to linear veloci- 
ties. We used the most r ecent determination o f the distance to 
LS 5039, 2.9 + 0.8 kpc (ICasares et al .11201 2al) . The distance to 
PSR J1825-1446 is 5.0 + 0.6 kpc. The first set of errors in the 



Galactic velocities in Table [7] corresponds to the proper motion 
uncertainty, and the second set of errors corresponds to the con- 
tribution of the distance uncertainty. 

We can now compute the Galactic space velocity as mea- 
sured from the Galactic Centre (GC). We need to transform 
the velocity vector to the reference frame of an observer sit- 
uated at the position of the source that participates in the 
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Fig. 5. Distribution of proper motion fits for LS 5039 computed 
with the bootstrap method. The central contours represent, from 
inside to outside, the minimum smoothed regions containing 
68.2, 95.4, and 99.7% of the events. The two histograms are the 
projected ID distributions shown in Fig. |4] 



Table 7. Galactic coordinates and Galactic space velocities of 
LS 5039, and PSR J 1825- 1446. In the Galactic velocities, the 
first set of errors corresponds to the proper motion uncertainty, 
while the second set (in parentheses) corresponds to the distance 
uncertainty. 



Parameter 


LS 5039 


PSR J1825— 1446 


Galactic Longitude / 
Galactic Latitude b 
Distance d (kpc) 


16°882 
-L289 
2.9 ±0.8 


16° 805 
-L001 
5.0 ±0.6 


fi a cos S (mas yr -1 ) 

He (masyr 1 ) 

Sys. velocity (y) (km s _1 ) 

fi\ (mas yr _1 ) 

fi b (mas yr 1 ) 


7 09+0.19 

o OO+027 
-8.8Z_ 012 

17.2 ± 0.5 
-10.38** 


10.0 ±0.3 
-29.0 ± 0.3 

-21.0 ±0.3 
-22.4 ± 0.3 


{/ L sr (km s _1 ) 
V L sr (kms- 1 ) 
W LSR (km s" 1 ) 


42.5 + | ° (±4) 
-43* (±17) 
-136+| (±40) 


146 ±2 (±16) 
-466 ± 7 (±60) 
-523 ± 7 (±60) 


Vrad (km S" 1 ) 

v cir (km s _1 ) 
v z (km s" 1 ) 
Vpcc (km s" 1 ) 


-11.2+1] (±5) 
197+3 (±17) 
-136+1 (±40) 
142+2 (±40) 


227 ± 3 (±40) 
-150 ±6 (±60) 
-523 ± 7 (±60) 

690 ± 7 (±60) 



Galactic rotation. This reference frame is the regional stan- 
dard of rest (RSR). The transformation from the LSR to the 
RSR takes into account th e distance from the Sun to the GC, 
8.0 + 0.5 kpc (Rei dlll993[). and the circul ar rotation about the 
GC, 236 km s 1 dReid & Brunfhaler|[2004l) . which we consider 
constant. From the Galactic velocity (U, V, W)rsr, we directly 



obtain the Galactocentric velocity of a source, with the radial 
velocity v ra d away from the GC, v c j r , in the direction of the 
Galactic rotation at each point, and v z , perpendicular to the 
Galactic Plane, and positive towards the NGP. The expected ve- 
locity of a system following the Galactic rotation is zero for v ra d 
and v z , whereas v c ; r is expected to be 236 km s _1 . The devia- 
tion from the expected circular rotation is the peculiar velocity 
Vpec - [vjL + ( v cir _ 236 km s~') 2 + v z ] I/2 . The Galactic veloci- 
ties from" the GC for LS 5039 and PSR J 1825 -1446 are shown 
in Table. [7] 

The large error in the distance of the sources limits the deter- 
mination of their space velocity. In this case it is convenient to 
plot the Galactic velocities as a function of the distance from the 
Sun. In Figs. [6] and [7] we show the computed velocities (black 
lines) and their uncertainties (dark grey areas) for a range of 
distances. The horizontal dashed lines mark the expected veloc- 
ity when following the Galactic rotation. The grey areas indi- 
cate the cosmic velocity dispersion of the Galactic young disc 
stars, which depends on the age of the stellar system. LS 5039 
is an O-type star and its companion could have been more mas- 
sive. The progenitor of PSR J1825-1446 is unknown, but col- 
lapse to form a neutron star us ually occurs for stars with to - 
tal mass between 8 and 15 M Q dLvne & Graham-Smithll2005l) . 
Therefore, we assumed that the progenitor was a young star 
formed in the Galactic thin disc. The measured velocity dis- 
persion of young OB stars is biased towards local values be- 
cause the stars with known 3D velocities and radial distances 
are predominantly located in the solar neighbourhood. However, 
the velocity dispersion of young stars should follow the veloc- 
ity dispersion of the gas in the Galactic disc. This dispersion 
is of about 4-5 km s in the plane (<tu,(Tv), and can be as 
low as 1-2 km s in the direction perpendicular to the plane 
(ctw), although spiral shocks can induce higher i n-plane veloc- 
ities of about ~ 1 — 1 km s' 1 ([ Kim et al. I2006h . This agrees 
with values measured in lTorra et al.l(l2000l) for stars of 0-30 Myr 
of {a-u,cr v ,a- w ) = (7 .9,7 .2, 4.3) km s _1 . We assume that the 
expected velocity dispersion of young stars is (cr u , cry, cr w ) = 
(10, 10,5) kms -1 . 

We can see in Fig. [6] that LS 5039 has a high peculiar ve- 
locity, mainly perpendicular to the Galactic plane, and it is not 
compatible with the expected velocity dispersion for any of the 
considered distances. This reaffirms the runaway nature o f the 
system, which was originally discussed in lRibo et al.l(l2002l) . On 
the other hand, PSR J1825-1446 shows a much higher space 
velocity. This is an isolated pulsar, and therefore it is not pos- 
sible to measure any radial velocity from the Sun. In Fig. |7] 
we show the Galactic velocities assuming a radial velocity of 
0, +500, and -500 km s _1 (black, red, and blue lines, respec- 
tively). This covers a generous range of possible radial veloci- 
ties. We note that this contribution is less notorious in the direc- 
tion perpendicular to the Galactic Plane because both the pulsar 
and the Sun are at a similar distance from the Galactic plane. 
For PSR J1825-1446 we plot a wide range of distances, con- 
sidering that the distance estimate could be overestimated by 
the DM determination. The pulsar space velocity is not com- 
patible with the expected rotation of the Galaxy. Remarkably, 
the pulsar has a velocity perpendicular to the Galactic Plane of 
~ 550 km s" 1 and away from it. Therefore, PSR J1825-1446 is 
clearly a runaway pulsar that has been ejected from the Galatic 
disc at a very high speed. The peculiar velocity for the nominal 
distance of 5.0 kpc is 690 ± 7 km s . This transverse velocity is 
high even for a runaway isolated pulsar, because it lays in the tail 
of the transverse ve locity distribution of young pulsars shown in 
Hobbs et al. (2005). The peculiar velocity would be lower if the 
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Fig. 6. Components of the Galactic and peculiar velocity of 
LS 5039 (solid black lines) and their uncertainties at 3-cr level 
(dark grey areas). The horizontal dashed lines and light grey ar- 
eas mark the expected Galactic velocity, and the velocity disper- 
sion expected for young stars, respectively. The vertical dotted 
lines indicate the measured distance to the source. 



LS 5039 and PSR J1825-1446 
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Fig. 7. Same as Fig. ©for PSR J 1825- 1446. The solid black line 
in each panel indicates the Galactic velocity assuming a radial 
velocity of km s _1 . The velocity uncertainties at 3-cr level are 
smaller than the black lines. The red thin line corresponds to an 
hypothetical radial velocity from the Sun of +500 km s (away 
from the Sun) and the blue thin line to -500 km s _I (towards the 
Sun). 



distance is overestimated, and would be 500 ± 5 km s 1 for a 
distance of 3.6 kpc. 

6. The origin of LS 5039 and PSR J1 825-1 446 

LS 5039 is a young 06.5 main sequence star, which has a total 
lifetime of a few million years. This sets an upper limit to the 
age of the binary system LS 5039/compact object, which allows 
us to constrain the possible locations where the supernova took 
place, and therefore where the system originated. We have ex- 
plored possible origins of LS 5039 in the last 1-10 Myr, even 
though for ages older than 1-2 Myr the source birthplace lies 
well outside the Galactic plane, which is unlikely. Therefore, we 
only consider the projected past trajectory of LS 5039 during the 
last ~ 10 6 yr, and discuss possible associations that cross this tra- 
jectory. We know that the characteristic age of PSR J 1825 -1446 
is 1.95xl0 5 yr. This value is computed as the ratio between the 
period and the period derivative of the pulsar, and it is valid un- 
der certain assumptions: the initial spin period of the pulsar was 
much shorter than that observed today (Pq <k P), there is no 



magnetic field decay (PP is constant), and the energy loss cor- 
responds to a spinning dipolar magnet in perfect vacuum. The 
slowdown model of a pulsar is generally parametrised by the 
braking index, which is usuall y considered close to 3, although i t 
can take a range of values (see lFaucher- Giguere & Kaspi 2006). 
The uncertainty on this parameter prevents the characteristic age 
from being a faithful estimate. It has been shown that the char- 
acteristic age of a pu lsar can be significantly lower than the true 
age of the pulsar (e.g. lGaensler & Fraill2000l) or ev en higher than 
the real age of the pulsar (e.g. lKramer et al.ll2.Q03b . Therefore, a 
wide range of ages, and consequently possible angular distances 
from the current location of the pulsar, should be explored. 

In Fig. [8] we show the region around the current position of 
LS 5039 and PSR J1825-1446 (green and red crosses, respec- 
tively) in Galactic coordinates. We also plot different objects re- 
lated with young stars. The ellipses in green, red, and blue mark 
the positions of the supernova remnants, the open clusters and 
the OB associations, respectively. To have a better context of 
the distribution of young systems, we include the O and B stars, 
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Fig. 8. Area around LS 5039 and PSR J1825-1446 in Galactic coordinates. The lines mark the past trajectory in the last 10 6 yr 
for LS 5039, and 1.95 x 10 5 yr for PSR J1825-1446. The green circles are supernova remnants from the catalogue compiled in 
lGreenl(l2009l) . T he red circles are open clus ters from the WEBDA database, and the blue ellipses m ark the position o f the known 
OB associations dMel'Nik & Efremovl 19951) . The OB stars (Reed 2003) and the star-forming regions (Avedisova 2002) are marked 
with blue and grey dots, respectiv ely, and provide a general view of the young system distribution in this region of the Galaxy. The 
red crosses mark the H n regions (lLockmanll 989). The objects in regions #1, #2, and #3 are described in the text. 



the star-forming regions, and the h n regions close in this part 
of the Galactic plane. The past trajectories in the last 10 6 yr for 
LS 5039, and 1.95xl0 5 yr for PSR J1825-1446 are indicated by 
the green and red continuous lines, respectively. The plotted tra- 
jectories are linear, and therefore not affected by the Galactic po- 
tential. We checked the validity of this assumption using the soft- 
ware package galp^\to compute the trajectories under a standard 
Milky Way potential (combination of a Miyamoto-Nagai disc, a 
Navarro-Frenk- White halo, and a Hemquist bulge potentials, see 
the documentation for details). For all aspects discussed in this 
paper the slight difference between the integrated orbit in the 
Galactic potential and the free (linear) trajectory are irrelevant. 
Therefore, we do not use the orbit integration to avoid using ad- 
ditional parameters to describe the Galaxy. We used the software 
to check that the space velocity of PSR J1825-1446 is above the 
escape velocity, and therefore it will be ejected from the Galaxy. 

As can be seen in Fig. [8] the past trajectory of LS 5039 
and PSR J1825-1446 crosses three interesting regions, marked 
with numbers from #1 to #3. The first region (#1) is com- 
posed of the objects at short angular distances from the sys- 
tems. In the WEBDA databasdj we can see that LS 5039 
and PSR J1825-1446 are within the open cluster Dolidze 28 



7 https : //github. com/j obo vy/galpy 

8 http://www.univie.ac.at/webda/ 



(C 182 2-146), whose central star is WR 115. lManchandaet al.1 
( 1996) determined the distance to the cluster to be 2.2 kpc (note 
that the cl uster is wrongly named Do 78 instead of Do 28 in 
Table 1 of iManchanda et alJll996h . We used the 2MASS photo- 
metric mag nitude of WR 115 (Wo lf-Rayet of type WN6) from 
Table Al in ICrowther et all (120061) . M Ks = -5.10, and colours 
for a nominal subtype WN6 star (J - K) = 0.17 and (H - K) = 
0.15. The absorption estimate of the source is Ak = 0.59, which 
yields a distance of 2.0 kpc, compatible with the distance quoted 
abov e. In the same field we fin d the OB association Set OB3 
(see iMel'Nik & Efremovl [1991 . with ten members in a region 
of 25.8x33.6 pc 2 , and an estimated dista nce of 1.5 kpc. Thi s 
distance is compatible with the estimate in Tetzl aff et al.l (|2010), 
obtained from a parallax measurement, of 1.4 kpc. They also 
provide an age of Set OB3 of 4-5 Myr. The closest H n ionised 
region is RCW 164, already discussed in Sect. 12. II 

In Fig. [9] we show the most relevant objects of re- 
gion #1 together with the past trajectories of LS 5039 and 
PSR J1825-1446. We can see that although PSR J1825-1446 
is currently very close to RCW 164, it comes from the opposite 
direction at high speed. The trajectory is not compatible with the 
central region of SNR G016. 8-01.1 either. Although it is not 
possible to determine the kinematic centre of the SNR precisely, 
the projected past trajectory of LS 5039 seems to be far away 
from the central region of the remnant (see Fig. [1}, although an 
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Fig. 9. Wide-field map around SNR G016.8-01.1. The grey 
scale corresponds to the Parkes-MIT-NRAO survey at 5 GHz 
dTasker et al.lll994h - The green and red lines mark the past tra- 
jectory of LS 5039 and PSR J1825-1446, respectively. The ar- 
rows mark the proper motion sense. Circles are plotted in the 
trajectories every 10 4 yr. The contours represent the 68.2, 95.4 
and 99.7% confidence levels on the position that LS 5039 had 
10 5 and 5 x 10 4 years ago (see uncertainties in the proper mo- 
tion shown in Fig. [5]). The grey dotted lines show the Galactic 
coordinates. 



association cannot be discarded. We note that the l-cr accuracy 
in the position of LS 5039 1 5 yr ago is approximately 18 arcsec, 
although the uncertainties do not follow a Gaussian distribution, 
as shown by the contours in Fig. [9] 

The past trajectory of LS 5039 shown in Fig. [9] suggests an 
interesting possibility, namely that the binary system LS 5039 
and WR 115 come from the same point on the sky, and that 
they might have been a triple system, from which WR 115 was 
ejected in the opposite direction after the supernova explosion. 
This would tightly restrict the age of the triple system. However, 
we have computed a proper motion of WR 115 using positions 
from optical catalogues, in particular from AC2000.2, GSC 1 .2, 
Hipparcos, UCAC3, and 2MASS, see Sect. I4AT1 The proper 
motion obtained is 4.5 + 0.6 and -5+4 mas yr in right ascen- 
sion and declination, respectively. The proper motion in right as- 
cension discards the possibility of a triple system with WR 115 
because the linear momentum conservation would require that 
the system pre-supernova already had a runaway velocity much 
higher than expected from Galactic diffusion. 

Region #2 is found at an angular distance of ~ 3° from 
both sources, and harbours an active site of star formation and 
young stars. The OB associa tion shown in region #2 is Ser OB2 
dMeFNik & Efremovl [l995l) . with a size of 22.8 and 33.0 pc 
in Galactic longitude and latitude, respectively, an average ra- 
dial velocity relative to the Su n of —2 . 4 + 3.9 km s , and an 
average distance of 1.45 kpc. Forbes) (120001) provides a dis- 
tance to the a ssociation of 1.9 + 0.3 kpc, and an average age 
of 5 + 1 Mvr. iTetzlaff et al] (2010) provide a parallax measure- 
ment based on Hipparcos stars of 0.63 mas, which correspond 
to a distance from the Sun of 1.6 kpc and sets an age of the 



association of 4-5 Myr. On the other hand, the open cluster in 
the centre of region #2 is NGC 6604, which can be found in 
the WEBDA catalogue, with a di s tance of 1.7 kpc, and an age 
of ~ 6 Myr. Loktin & Beshenov (2003) measured an average 
proper motion of this cluster, formed by ~ 97 stars, and obtained 
-0.66 + 0.19 and -0.88 + 0.17 mas yr 1 in right ascension and 
declination, respectively, as well as the distance quoted in the 
WEBDA database. Recently, Tetz laff et al.l d2010t) set a parallax 
of 0.59 mas, corresponding to a distance from the Sun of 1 .7 kpc, 
and an age of 4-5 Myr. There is also an interesting H n region 
that might be associated wit h this complex, SH 2-54, from the 
catalogue of H n regions in Sharp lessl (119591) . Its radial velocity 
was measured bv lBlitz et al.l (ll982) to be 27.6±0.5 km s , much 
higher than the average velocity of Ser OB2 quoted above, and 
the distance from the Sun is 2.0 + 0.2 kpc, which is compatible 
within errors with the structures discussed above. In summary, 
region #2, which harbours many young stars and star-forming 
regions, is undergoing intense formation of new stars, it is at 
a distance of 1.5-2.0 kpc from the Sun, and has an age of 4— 
5 Myr. The age of the complex is enough to allow LS 5039 to 
have travelled from that region of the sky to its current position. 
Considering the large uncertainties in the distance to the system, 
this is a plausible birth region for LS 5039 provided that it is 
located between 1.5 and 2.0 kpc from the Sun. For this distance 
range, Fig. [6] shows that LS 5039 would have a lower peculiar 
Galactic velocity of about 70-95 km s _1 . In this case, LS 5039 
would have an age between 1 .0 and 1 .2 Myr. 

Finally, region #3 in Fig. [8] marks the position of 
SNR G017.4-00.1, which has a diame ter of 6', a flux density of 
0.4 Jy , and a spectral index of +0.7 dBrogan et al. 2 0061 iGreenl 
2009). This partial shell remnant is of class II, which corre- 
sponds to a gravitational collapse supernova explosion. The min- 
imum projected distance between the past trajectory of LS 5039 
and this SNR took place 4 x 10 5 yr ago, and for PSR J 1825 -1446 
it took place 0.13 X 10 5 yr ago. For these high ages we do not 
expect to still detect the SNR, consequently it is unlikely that 
PSR J 1825 -1446 and SNR GO 17. 4-00.1 are related, whereas it 
is not possible tha t LS 5039 and S NR G017.4-00.1 are related. 
On the other hand, Bochowl (1201 ll) used the 2 - D relationship to 
estimate the distance to SNR G017.4-00. 1, and obtained 12 kpc. 
The distance to the SNR seems considerably larger than the dis- 
tance to the sources, although we note that this estimate can be 
very uncertain. 

Besides the specific regions on the sky in the past trajec- 
tory of both systems, we can constrain the lifetime of the sys- 
tems considering that the systems had contained young mas- 
sive stars, which are generally formed in a specific region of the 
Galactic plane, the thin disc. The stellar space density decays 
exponentially as a function of the distance above or below the 
Galactic plane, and each stellar population h as a particul ar scale 
height h. For young O B stars /zob = 45 pc (Reed 2000), while 
Mafz-Apellani^ (1200 ll) found that the 0-B5 stellar population is 
distributed with a scale height of 34.2 + 0.8 + 2.5 pc. No good 
determination for the O-type population alone has been found, 
basically due to the lack of good statistics with this population, 
and therefore we used the more conservative value of 45 pc both 
for LS 5039 and PSR J1825-1446. From the Galactic latitudes 
quoted in Table|7] the current height with respect to the Galactic 
plane is -65 and -96 pc forLS 5039 and PSR J 1825 -1446, re- 
spectively. With the proper motion of the systems, we compute 
for different distances from the Sun the time when they were 
crossing the Galactic plane at a height within ±45 pc. The ex- 
pected ages are quoted in Table [8] 
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Table 8. Crossing times of LS 5039 and PSR J1825-1446 
through the Galactic thin disc, with a scale height for OB stars 
of +45 pc, for different distances from the Sun. 



LS 5039 


PSR J 1825- 


1446 


Distance [kpc] Age [kyr] 


Distance [kpc] 


Age [kyr] 


2.0 0-890 


3.0 


20-300 


2.9 140-750 


5.0 


80-245 


4.0 220-670 


7.0 


100-220 



7. Discussion and conclusions 

We have obtained a set of accurate positions of the isolated pul- 
sar PSR J1825-1446 by means of VLBI observations. We used 
the pulsar gating technique to enhance the detected pulsar flux 
density to improve the astrometry. We fitted a linear proper mo- 
tion to new and archival astrometry and obtained a proper motion 
of n a cos 5 = 10.0 + 0.3 mas yr , and /ug = -29.0 + 0.3 mas yr -1 . 
For a distance from the Sun of 5.0 kpc, the proper motion of 
PSR J1825-1446 corresponds to a projected transverse velocity 
of 690 + 7 km _1 , which makes it a high-velocity pulsar ejected 
from the Galactic plane. 

We note that the high space velocity of PSR J1825-1446 
suggests that it is possible that a bow-shock has formed in front 
of the pulsar. For example, the pulsar wind nebula G359.230.82 
is the bow shock produced by the young pulsar J1747-2958, 
which has a spin-down luminosity of 2.5 x 10 36 erg s , and 
a transverse velocity of 300 km s~_ in a medium w ith hydro- 
gen number density of 1.0 cm~ 3 dHales et al.l [20091) . However, 
the spin-down luminosity of PSR J1825-1446 is only 4.1 x 
10 34 erg s _1 , and therefore a such structure would only be visible 
if the medium surrounding PSR J 1825 -1446 were much denser 
than in the case of J1747-2958. 

The space velocity of the source has allowed us to inves- 
tigate possible origins of PSR J1825-1446. The characteristic 
age of PSR J1825-1446 is approximately 0.2 Myr. The ve- 
locity and age of PSR J 1825- 1446 make it incompatible with 
SNR G016.8-01.1. There are no clear OB associations or SNRs 
crossing the past trajectory of the pulsar. Our estimate of the age 
of the pulsar, assuming that its progenitor died in the Galactic 
disc, is 80-245 kyr, compatible with its characteristic age. We 
have seen that the space velocity of this pulsar is high enough to 
escape from the Galactic potential. 

We compiled all radio interferometric observations of 
LS 5039 with accurate astrometry, and reduced all data in a 
consistent way. Combining the radio positions with archival op- 
tical astrometry from global catalogues we obtained a proper 
motion of the source of fi a cos6 = 7.09^'jj' mas yr -1 , and 
/jg - -8.82^q j2 mas yr -1 . We have reduced the uncertainty on 
the proper motion by approxi mately one order of magnitude with 
respect to the one obtained in Rib o et al.l (2002), which is com- 
patible at 2-cr level. The data were collected from very hetero- 
geneous observations (completely different instruments, condi- 
tions, frequencies, and reference calibrators), and the fit is af- 
fected by systematic errors. We investigated the effects of these 
unknown uncertainties using a bootstrap method, which pro- 
vides realistic uncertainties of our measurement. We confirmed 
the runaway nature of the system, which has a peculiar space ve- 
locity of 141 ^ km s _1 , with amain component of -136^ 3 km s _1 
perpendicular to the Galactic plane, for a distance to the source 
of 2.9 kpc. 



We we unable to identify a secure origin for LS 5039. We 
restricted the proper motion of the system, and the past trajec- 
tory is not far from the putative centre of SNR G016. 8-01.1, 
which is very uncertain. However, this association would im- 
ply an age of 10 5 yr for the SN R, which is above the expected 
lifetime for such a bright SNR dFrail et al.l 11994b . The age of 
the SNR depends on the ambi ent density and it s distance, which 
are unknown parameters (see iRibo et aLll2002l) . However, even 
if SNR G016.8-01.1 is not related to LS 5039, the entire re- 
gion shown in Fig. [9] is an active region of stellar formation, 
with the open cluster Dolidze 28, WR 115, and the Hn region 
RCW 164. It is possible that the compact object was formed 
within this complex. In this case, LS 5039 should be at a dis- 
tance of around 2 kpc, and the age of the binary would be be- 
low 0.1-0.2 Myr. On the other hand, the region of Ser OB2 lies 
within the past trajectory of LS 5039, at a distance from the Sun 
of 1 .5-2.0 kpc, which is still compatible at 2-cr level with the dis- 
tance to LS 5039, 2.9+0.8 kpc. In this case the age of the system 
would be 1.0-1.2 Myr, and the peculiar Galactic velocity would 
be 70-95 km s" 1 . A third possibility is that the system LS 5039 
was formed in the Galactic plane, but it was not associated with 
any known young stellar association. In this case the age of the 
system would be between 0.1 and 0.8 Myr, assuming a distance 
of 2.9 kpc. 

The age of LS 5039 has consequences on the nature of 
the compact object. When the compact object was formed, 
the system received a powerful kick, and therefore the orbit 
could not be pseudo-synchronised after the explosion. Pseudo- 
synchronous rotation is reached through tidal forces in an ec- 
centric binary, when the rotation of the massive star and the mo- 
tion of the compact object at periastron are synchronised (lHalll 
1986). Casares et al.l (120051) determined that the compact object 
should be a black hole (low system inclination and high compact 
object mass) if the orbit of the system is pseudo-synchronised. 
They estimated the time for this system to reach orbital pseudo- 
synchronism to be ~ 1 Myr, although it is highly uncertain. If 
the system was formed in Set OB3 region, its age would be be- 
low 0.2 Myr, and the hypothesis in ICasares et al. 1 (12001 would 
be unlikely. Only if the source was formed close to Ser OB2 re- 
gion, 1-2 Myr ago, would it be possible to reach orbital pseudo- 
synchronism. Based on the proper motion of the source, it is not 
expected that the system age is greater than 1-2 Myr if the sys- 
tem was born in the Galactic disc. 

We have seen that any determination of the age and ori- 
gin of LS 5039 is limited by the uncertainty in its distance to 
the Sun. This distance cannot be obtained by means of paral- 
lax VLBI measurements with the current resolution, and the ef- 
fects of the source structure and variability . An accurate dis tance 
will be obtained by the Gaia satellite (see Ide Bruiinal2012l and 
references therein) through measurements of the parallax of the 
massive companion. The satellite will also determine the dis- 
tance to the OB associations with high accuracy, and therefore 
any identification or association should be much more straight- 
forward, and will help to clarify the origin of the system. On 
the other hand, the distance to PSR J1825-1446 could be ob- 
tained by means of a challenging long-term VLBI project. This 
would provide a model-independent distance that can be com- 
pared to the estimated one from the dispersion measure. Finally, 
there is no obvious approach to accurately determine the dis- 
tance to SNR G016. 8-01.1, whose radio structure is severely 
contaminated by the H n region RCW 164. 
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to the regular pulsar monitoring conducted in the Jodrell Bank 
Observatory. Two sets of data were obtained for each observa- 
tion, one from a normal correlation and one using pulsar gating. 
This allowed us to compare the results obtained with and with- 
out the pulsar gating correlation. The three VLBA observations 
were correlated with a gate covering pulsar phases from 0.99056 
to 0.01389, integrating 6.5 ms of each pulsar cycle. 

The same calibration tables applied to the normal data were 
also applied to the data correlated with pulsar gating. We pro- 
duced images of PSR J 1825 -1446 using the same parameters, 
and we compared the astrometry obtained with and without 
pulsar gating. The obtained signal-to-noise ratio and the cor- 
responding gain factor improvement are shown in Table [3] On 
average, the obtained gain of 2 is relatively small. This can be 
explained by an asymmetric pulse profile, or an inaccurate de- 
termination of the pulsar gate. 



Appendix A: Pulsar gating on PSR J1 825-1 446 

To improve the sensitivity of the detections, and thus the as- 
trometry, we correlated the three VLBA observations using pul- 
sar gating (Bri sken et al.l2002T ; IChatteriee et al.l200 9). In normal 
VLBI observations, with integration times of a few seconds, one 
measures the period-averaged flux density of the observed pul- 
sar, although most of the time the pulsar is not emitting. The 
instantaneous flux density during the pulsation is usually tens of 
times higher than the period-averaged value, because the duty 
cycle (the ratio of the pulse width to the pulse period) is usually 
below 10%. To optimise radio observations of pulsars one can 
take advantage of this fact by only correlating the data during 
on-pulse, and disabling correlation during off-pulse. The total 
correlated on-source time is considerably reduced (a few min- 
utes for every observed hour), which increases the final image 
noise (rms) by a factor proportional to one over the square root 
of the duty cycle (At v ) at a certain frequency v. However, the 
flux (S v ) increase is approximately proportional to one over the 
duty cycle, yielding a final signal-to-noise ratio (S /N) gain: 

(S/AOc = I*L ~ = (S/N)' (A.1) 

rms G (rms)/ VAtv y[Kr v 

where the subindex G stands for the gated values. This pro- 
vides a theoretical signal-to-noise ratio increase proportional 
to one over the square root of the duty cycle of the pulsar. 
Typically, the real signal-to-n oise improves by a factor between 
2 and 6 (Bris ken et al.l 120021) . which can make the difference 
between a non-detection and accurate astrometry. The duty cy- 
cle of PSR J1825-1446 is 2.5% (6.8 ms) at 4.8 GHz, yield- 
ing a maximum theoretical signal-to-noise gain of 6.3. For the 
three VLBA correlations, updated ephemeris of the pulse time 
of arrivals were obtained to configure the opening and closing 
time "gates" of the correlator. The timing was obtained thanks 
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